Moments of the ARPES spectral function of an undoped Mott insulator. 
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We derive analytic expressions for the first three frequency moments of the single particle spectral 
function for one hole in a Mott insulator in terms of equilibrium spin correlation functions of the 
insulating state. We show that the "remnant Fermi surface" detected in ARPES experiments is, in 
fact, a reflection of the strong antiferromagnetic correlations of the system, not a reflection of the 
original band-structure Fermi surface. We suggest that ARPES data could be used to measure the 
magnetic specific heat. 



In this paper we analyze the single hole spectral func- 
tion of a Mott insulator. We take as our model of the 
insulator the large U Hubbard model at half-filling and 
find that spectral moments can be evaluated using stan- 
dard perturbation theory in the kinetic energy. Equilib- 
rium magnetic correlations of the half-filled state fully 
characterize the one hole spectral moments. 

Even well understood strongly correlated electron sys- 
tems have non-trivial single particle properties. Loosely 
speaking, when the effects of interactions, kinematics, 
and/or magnetic field overwhelm the kinetic energy of 
the electrons, the resultant new states of matter bear 
little or no resemblance to the non-interacting electron 
gas. And while sharp single particle spectral peaks are 
the defining property of weakly correlated systems (met- 
als, semiconductors, even BCS superconductors) broad 
incoherent features, often alongside a weak quasi-particle- 
like (QP) peak, are more common in strongly correlated 
materials. Specifically, angular resolved photoemission 
spectroscopy (ARPES) data taken from Ca2Cu02Ci2 (a 
prototype Mott insulator) does showEJ a peak in a lim- 
ited region of phase space but is clearly dominated by 
broad features. Numerical solutions on small systems 
combined with judiciously summed perturbation schemes 
successfully reproduce some aspects of photoemission 
resultso ! Q0. However the underlying physics determin- 
ing the spectral function is often obscured; even more of- 
ten the approximation scheme is not universally accepted 
as justifiable. In contrast, by asking different (though 
related) questions we have not only drawn a clear con- 
nection between one hole properties with those of the 
undoped antiferromagnet (AF), but have done so in a 
controlled manner. The basis of our calculation will be 
reviewed in Section 2. In Section 3 we will present and 
discuss our result for the zeroth moment, n^. In par- 
ticular, we will show that "a remnant Fermi surface," 
defined as either the locus of points in k space at which 
= 1/2 or as the locus of points of maximum gra- 
dient of 7i-, lies close to the "diamond Fermi surface" 

k 

(cos(k x ) +cos(fcj,) = 0), regardless of whether or not this 
is the location of the Fermi surface in the U — > limit. 
Results concerning higher moments are presented in Sec- 
tion 4. Some of the technical details can be found in the 
Appendix. 



I. THE CALCULATION 

A. Half- filled Hubbard model at large U 

The large U Hubbard model at half filling is the sim- 
plest model that produces the large charge gap and the 
antiferromagnetic tendencies observed experimentally in 
CuO planes. The Hamiltonian, 
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contains a strong onsite repulsion term U ^P^nZn- a that 

3 3 

freezes (together with the lattice potential) charge mo- 
tion at half-filling provided that the kinetic energy term 
is small enough to be treated as a perturbation, |ty | <C U. 

When tij/U — the ground state consists of singly oc- 
cupied sites, has energy zero, is 2 7V -fold degenerate where 
N is the size of the system, and is completely character- 
ized by the spin configuration, {Si}. 

The perturbative effects of the kinetic energy modify 
all of the above statements except the last one: each of 
the perturbed states, though it contains an admixture 
of doubly occupied sites, can still be labeled by the spin 
configuration of the unperturbed state from which it has 
evolved. Thus the expectation value of any operator in 
the ground state manifold can always be expressed in 
terms of spin variables. Formally this is accomplished by 
computing perturbatively in powers of tij/U the unitary 
transformation exp[iJf] that expresses the evolution of 
the low energy unpeturbed..states \{Si} > as a function 
of increasing kinetic energyQ: 



\{S i }> = e- iJt \{S i }> 
< 6 > = < {Si}\e iA 6e 



>x 



\{Si} >, 



(2) 
(3) 



where O is any observable. 

In the familiar fashion, this transformation maps the 
low energy physics of the Hubbard Hamiltonian into an 
effective Heisenberg antiferromagnet whose leading term 
is: 



H, 



eff 
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Si ■ S 3 ) + 0(4/C/ 3 ). (4) 
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B. The spectral function 

The emission spectral function is 
A(k, w) = l-Y e-^" | < m\C\ \n > \ 2 6(iu + E m - E n ) 

<7,m,n 

= ^J2j dteiwt<c l^ c nJ°)>- ( 5 ) 

Frequency moments of A(k, w) correspond to ground 
state (or thermodynamic) averages of the following oper- 
ators: 



Ax(k) 



dw 
2^ 



u A{k, w) 



A 2 {k) = 



dw 
2^ 



iv 2 A(k,w) 



(7) 



(8) 



A systematic evaluation of these averages, using the 
perturbative expression for X in powers of Uj/U , can be 
found in the Appendix. 



II. THE SPECTRAL WEIGHT, n g 

We begin with a discussion of the spectral weight be- 
cause it has an immediate physical interpretation as the 
occupation probability. As shown in the Appendix, 
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where 



and 



Si 



ik-Rj 



- S t ■ S, > 
4 J 



(9) 



(10) 



(11) 



is the equilibrium spin correlation between spins i and j. 
ik is a sort of renormalizcd band energy in which each 
hopping matrix element is renormalized by a factor of 
Sij. However, this renormalized energy does not corre- 
spond in any simple way to the energy of any elementary 
excitation of the system. Note that the non-interacting 

free electron band is given by = — J2j toje lh ' Rj ■ 

This expression for can be derived in a different, 
simpler manner, which is readily generalizible to more 



complicated situations, such as the three band Cu-0 or 
Emery modeB. From the Hcllman-Feynman theorem, it 
follows that 



E 



< Wi« + H - c -l >= -dE/dta = [dJ %3 ldt l3 ]S l3 



(12) 



where E is the internal energy, which can be computed 
using the effective Hamiltonian H e ff. From the expres- 
sion for Jij in terms of tij, the result in Eq.|9| follows 
immediately. 



A. The remnant fermi surface 

At low temperatures, the short-range spin correlation 
functions are essentially temperature independent, and 
equal to their value in the ground state. To be con- 
crete, let us consider the Hubbard model with nearest, 
second, and third neighbor hopping, t, t' and t" , respec- 
tively; this sort of model was used in the numerical stud- 
ies to fit the dispersions seen in ARPES. The zero tem- 
perature spin correlations of the corresponding spin 1 /2 
Heisenberg model have been computed fairly accurately 
in numerical studiesQ. For the Heisenberg model with 
only nearest-neighbor exchange coupling, the spin cor- 
relation functions are Sot ~ 7/12 and S03 « £02 ~ 
between nearest, next-nearest, and third nearest neigh- 
bor sites, respectively. These correlations are, more- 
over, found to be relatively insenstitive to the inclusion 
of a modest amount of further neighbor exchange cou- 
plings, which anyway are expected to be quite small since 
J'/J=[t'/t} 2 . 

In computing ek the antifcrromagnetic correlations be- 
tween neighboring spins imply a factor of 1/2 renormal- 
ization of t, compared to a factor of 1/20 renormaliza- 
tion of t' and t". Since in most cases of physical interest, 
\t'\, \t"\ <C 10|t|, even when t' and t" are large enough to 
make signicant shifts in the original Fermi surface defined 
by efe, the occupation probability is well approximated 
(Figure 1) as 



1 7t(cos(k x ) + cos(k y )) 
2 [ 3U J 



(13) 



It is both i and has the steepest slope along the Fermi 
surface of the non-interacting electrons with only the 
nearest neighbor hopping. In fact, ARPES experiments 
that observe such momentum dependence of n- have 
been interpreted as an indication that there is a "rem- 
nant Fermi surface" in the undoped Mott insulator. By 
contrast, our result suggests that the observed n- is re- 
flective of the the spin physics of the strongly correlated 
Neel state rather than a vestige of the original Fermi sur- 
face. 
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FIG. 1. rir with t' = —.3t,t" = .15t and spin corre- 
lations Si = 7/12 and S 2 = Sa = l/20.Thick line is 
the non-interacting Fermi surface, while thin line is where 

"S = V2 

B. Specific heat of the Neel transition 

The natural connection we find between the spectral 
weight and the Neel state can be exploited further. Since, 
as we already pointed out, the Heisenberg Hamiltonian is 
dominated by the nearest neighbor term, and assuming 
that J n . n . is only weakly (if at all) temperature depen- 
dent, the specific heat is 

Ft 1 B f 

C(T) = J-< l -S 1 .S >=J^. (14) 

Since 

1 g^cos^) +cos(fc y )) 

n fe~2^ e7 lJ ' ^ J 

by measuring the temperature dependence of Si (as ex- 
tracted from n~) and differentiating it one gets the mag- 
netic contribution to the specific heat. The specific heat 
extracted in this manner contains only the contribution 
from the spin fluctuations. 

III. HIGHER MOMENTS 

Higher spectral moments correspond to the derivatives 
< ct {t)Cj; (0) > at t = and thus provide further in- 
sight into the problem of a hole in a Mott insulator. Our 
results for first and second moments, A\(k) and A 2 (k), 
are presented in a table below. The moments are written 
as sums over momentum space Fourier harmonics, each 



corresponding to a summation over sites a given near 
neighbor distance away on the square lattice. 

A n {k) = Y J A r } l3 (k) 1 

3 

iM= E ( 16 ) 

R—j th n.n. 
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A i 


^2 


1 


f(s-Si) 

2 f'(S 2 +2S 1 ) 
_«1(S 3 + 2S 1 ) 


-(2tf + tt")(Si - 1) 


2 


-t'{\-s 2 ) 

^"(S 3 + 2S 2 ) 
~^(S 2 + 2Si) 


-(2t't" + t')(S 2 -l) 


3 


-*"(|- so 

-£f(S 3 + 2S 2 ) 
-l n (S 3 + 2S 1 ) 


-((tr+^-i) 


4 


— \[-{Sa + S 2 + Si) 

— ^jj-(Sa + S3 + Si) 


-(t't" + tt"){S 4 -l) 


5 


-ini(s 5 + 2s 2 ) 

-^(S 5 + 2S 3 ) 


-(2£ + (t") 2 )(s 5 -i) 


6 


-t-fiSe + Ss + Si) 


-tt"(S e - 1) 


7 


-±f(S 7 + S 3 + S 2 ) 


-t't"{S 7 - 1) 


10 


-^rf(S w + 2S 3 ) 


_^ (Sl0 _i) 



One notices that all computed moments are finite. Al- 
though the existance of the moment expansion is a gen- 
eral requirement of any physical system, it is more rule 
than exception that approximations lead to divergencies 
past some finite order. Though we haven't constructed 
an explicit proof, there are indications that the expansion 
is well behaved for the Hubbard model. 

In principle, the moment expansion can be used to 
study the QP dispersion (Eg) directly: a coherent oscil- 
lation results in E'S contribution to A n (k). The crud- 
est (single mode) approximation of this sort identifies 
Eg = A\(k)/ng. We found it to be in a surprising agree- 
ment with previously obtained results for the momentum 
dependence of the QP energy in t-t'-t"- J as well as t- 
J models (one needs to assume that term proportional 
to t does not contribute to QP dispersion). Since Eg 
rather seriously overestimates the overall bandwidth, it 
isn't clear if one is justified claiming to have obtained 
even an approximate QP energy yet. Another likely use 
of our results (or rather their extension to higher mo- 
ments and orders in perturbation theory) can be in com- 
paring with spectral functions obtained by other means 
(cither numerics, self-consistent Born approximation or 
other). 

In conclusion, we have outlined a well controlled 
method for analysing the spectral moments of a hole in a 
Mott insulator. We find the occupation probability, n-, 
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is in agreement with the well established experimental 
result, which as our calculation suggests is strongly con- 
strained by the presence of AF order. We further propose 
that the temperature dependence of n- can be used to 

k 

study the specific heat of the Neel transition. The impli- 
cations of our results for higher moments are yet to be 
understood properly. 
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APPENDIX A: EVALUATION OF MOMENTS 

To compute the moments we first express them in 
terms of real space electronic correlations 

a 

= l^ + J2< C l, C ^>7i(k)} (Al) 

tr 

= tkn k + - 7i(*) < C L C o,* n v,-° > ( A2 ) 

2,t7 

Mk) = l^<CtjH,[H i C j ;J> 

U 2 

= ehn % E % + — < d^no.-a > 

+ 7i(*)*i,o < C { a (- C l* n o,-° 

+ C 1,-A,-<, C S,° + C) _C^_ a C^ > (A3) 

As before 7»(fe) = J2r e~ lk ' Ri , where the sum is over i'th 
nearest neighbors. 

Next, these electronic correlations are evaluated per- 
turbatively {S lj =< \ - Si ■ Sj >): 

E < C i C 5,« > = ^ + °(' 3 /t/ 3 ) (A4) 

(7 



On a square lattice and in a state that doesn't break 
time reflection invarianceQ (< §i ■ {S Rj x So) >= 0) the 
second sum is (for different i): 



1 


Att'S 2a + 2tt"S 30 


2 


2(251 - S 2 )t 2 + 4S 3 t't" 


3 


(2S 1 - S 3 )P + 2(2S 2 - Sz)t rz 


4 


2(5i + S 2 - Si)W + 2(5i + ^3 - Si)tt" 


5 


{2S 2 - S 5 )t'* + 2(2S 3 - S 5 )t"* 


6 


2(5i + S 3 - S 6 )tt" 


7 


2{S 2 + S 3 - S 7 )t't" 


10 


(2S 3 - S 10 )t'" 



Substituting these terms into A1,A2,A3 yields results 
of Section III. 
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